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Abstract: A research prototype CT scanner is currently under development in our lab. One 
of the key components in this project is the CT detector. This paper describes the design 
and performance evaluation of the modular CT detector unit for our proposed scanner. It 
consists of a Photodiode Array Assembly which captures irradiating X-ray photons and 
converts the energy into electrical current, and a mini Data Acquisition System which 
performs current integration and converts the analog signal into digital samples. The 
detector unit can be easily tiled together to form a CT detector. Experiments were 
conducted to characterize the detector performance both at the single unit level and system 
level. The noise level, linearity and uniformity of the proposed detector unit were reported 
and initial imaging studies were also presented which demonstrated the potential 
application of the proposed detector unit in actual CT scanners. 

Keywords: computed tomography; X-ray detector; performance evaluation; uniformity 
correction 



1. Introduction 

X-ray computed tomography (CT) is a medical imaging technology that utilizes computer-processed 
X-ray projections to produce tomographic images or "slices" of specific areas of the body. These 
images are used for diagnostic and therapeutic purposes in various medical applications. CT was 
invented in 1972 by Godfrey Hounsfield and the first clinical CT images were produced at the 
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Atkinson Morley Hospital in London in 1972 [1]. The original systems were dedicated to head 
imaging [1], and "whole body" systems with larger patient openings became available in 1976 [1]. CT 
imaging is widely used in the diagnosis of head, lung, cardiac, abdominal and pelvic diseases. An 
estimated of 72 million scans were performed in the United States in 2007 [2]. 

Because of the high system complexity and material cost, CT scanners are only available from a 
few manufacturers including Siemens, GE, Philips and Toshiba. Despite the development of CT 
technology, there are very limited choices in terms of CT components and instruments for the research 
community. In our lab, we proposed to build a modular CT scanner for the research purposes. One of 
the key components in this project is the CT detector. This paper describes the design and development 
of the modular CT detector unit for our proposed scanner. Evaluations were conducted to characterize 
the detector performance and results obtained in these experiments are reported. Initial imaging studies 
were also presented to demonstrate the functionality of the modular CT detector unit. 

2. Materials and Methods 



2.1. Modular Detector Unit 



Figure 1 shows the functional diagram of the proposed modular CT detector unit. It consists of a 
Photodiode Array Assembly (PDA) and a mini Data Acquisition System (mDAS). The PDA captures 
irradiating X-ray photons and converts the energy into electrical current. The mDAS mainly performs 
current integration and converts the analog signal into digital samples and output to follow-on circuits 
for further processing and data transfer. In our design, the ADAS 1128 from Analog Device (ADI) is 
used to facilitate the development of mDAS, which is a highly integrated 128-channel, simultaneous 
sampling and current-to-digital converter ASIC [3]. The mDAS is directly attached to the back of the 
PDA through a high density connector. Figure 2 shows a photo of the modular detector unit. 

Figure 1. Illustration diagram of the detector unit. 
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Figure 2. Photo of the modular detector unit. 




Table 1 lists the main characteristics of the PDA used in our design. The PDA consists of a stacked 
scintillator array and photodiode array. The scintillator material is gadolinium oxysulfide (Gd202S) 
with a peak wavelength of 512 nm. The scintillator emits light photons when irradiated by X-ray 
photons. These light photons are then converted to electrical charges in the readout photodiode array. 
The photodiode array size is organized as 24 x 24 (channels x segments) with 16 small pixels in the 
middle and four large pixels on each side along the segment direction. The pitch size is 0.95 mm and 
1.90 mm for the small pixels and large pixels, respectively. 



Table 1. Properties of photodiode array assembly. 


Property 


Condition 


Value 


Photosensitivity (mA/W) 


@510 nm 


320 (typical) 


Dark current 


small pixel 


10 pA 


large pixel 


20 pA 


Array size 




24x24 


Pixel pitch 




0.95 mm/1.9 mm 


Scintillator material 




Gd202S 



The proposed detector unit supports variable pixel size with several readout modes and Figure 3 
shows four of such modes. As shown in the figure, the detector outputs the 16 small pixel segments in 
the middle in mode 0, while in other mode, specific photodiode pixels are combined together as one 
readout segment. For example, in mode 1, every two small pixels in the middle are combined as one 
large pixel, together with the eight large pixels on the two sides, the detector outputs 16 large pixel 
segments. The detector also supports a float mode which disconnects the PDA from the mDAS, which 
is intended for test purposes. It is worthwhile to note that although the PDA is organized as a 24 x 24 
array, a maximum of 16 segment outputs along the channels are supported so that the detector unit has 
a maximum of 384 channel outputs, such as in the mode 0 and mode 1. 



Figure 3. Detector unit readout modes. 
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2.2. Experiment Setup 
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The performance of the proposed detector unit was characterized at both the unit level and the 
system level. All the experiments were conducted in a lead-shielded and light tight room. Figure 4 
shows the test setup using a single detector unit. The detector unit was placed in an aluminum box 
which acts as the shield for the detector unit. The lid of the aluminum box was removed for minimum 
attenuation of the irradiating X-ray photons. In order to make a dark environment for the detector unit, 
the aluminum box opening was sealed with black light-tight tapes. The aluminum box was mounted on 
a moving stage for vertical alignment of the detector unit to the focal spot of the X-ray tube (model 
G297 from Varian). An Indico 100 X-ray generator from CPI was used to control the generation of 
X-ray photons during the tests. In order to read out data from the detector unit, a custom designed 
circuit board was made which converts the data output from the detector unit into RS232 protocol and 
transfers the data to a host PC, and a MATLAB GUI tool was developed on the host PC to receive data 
from the detector unit and store as files for further analysis. 



Figure 4. Photo of the modular detector unit. 




Figure 5 shows the test setup of a prototype CT detector system. The detector system consists of 40 
detector units mounted side by side along an arc rail. The detector system was installed on an optical 
table with the arc center points to the focal spot of an X-ray tube. A high speed rotation stage was 
placed between X-ray tube and detector, which was used to mount imaging objects. The setup was 
used to simulate the gantry rotation in an actual CT scanner. A custom designed data acquisition 
control and management system was developed to read out the data from the CT detector and transfer 
the data to a host PC through a 1.25 Gbps optical fiber. The host PC installs a custom made data 
communication board to receive data from the detector and save them to hard disk as raw data. The 
raw data were then reconstructed into image slices using fan-beam filter-back projection (FBP) method. 
No data correction were performed during the tests reported in this paper. 
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Figure 5. Setup of the CT detector tests. 
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In the following text, experiments carried out at the unit level were denoted as testl and 
experiments carried out at the system level were denoted as test2 respectively. Detector unit's readout 
mode is set to mode 0 or float mode both in testl and test2 setup. 

2.3. Detector Performance Characterization Methods 

2.3.1. Noise Level 

There are three sources of noises in the proposed detector unit: PDA noise, denoted as n^ , which 
includes dark current from the photodiodes and thermal and electronic noises from the PDA circuitry, 
the circuit noise of the mDAS, denoted as n^, and the sampling noise of the ADAS 1128, denoted as 
n^ . Assume that n^ , n^ and n^ are independent, the measured noise of the detector unit, denoted as n , 
can be written as: 



(J, =(7^ +<J^ +(J^ 



(1) 



where a stands for the variance of the measurement. The sampling noise of the ADAS 11 28, n^ , can 
be obtained from the manufacturer. When set to the float mode, the PDA will be disconnected from the 
detector, in which case n^ is excluded and the measurement will only include n^ and n^ . So we can 
obtain two sets of measurements, one with n^ and one without n^ . Using these two measurement sets, 
n^ and n^ can be evaluated independently. 

In readout mode 0 and float mode, 384 analog signals are sampled simultaneously from a single 
detector unit and the acquired N frames of data can be denoted as a 384-dimension vector series 
X{n) = [xQ(n),x^(n),L ,x^^^(n)f . Noise level of the / 's channel is evaluated by its standard deviation, 
which defined as follows: 



(2) 
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When is large enough, (j. is an excellent estimation of noise. Power spectrum of measurement is 

estimated through the Periodogram method. Actually, the Periodogram is the square of signal's DFT 
which defined as follows [4]: 

2 



N-\ 



n=0 



Ik 

, fr = exp(-7 — ) (3) 



The experiments were carried out using a single detector unit with the setup described in Section 2.2. 
96 frames of data were recorded with the detector unit set to float mode and mode 0 respectively. An 
integration time of 500 us was used and the X-ray was turned off in these tests. Table 2 lists the 
conditions for the noise tests. 

Table 2. Noise level experimental setup. 



Readout mode 


Integration time T (us) 


Frames A'^ 


X-ray 


float 


500 


96 


off 


mode 0 (small pixel) 


500 


96 


off 



2.3.2. Independence 

Independence refers to independence of different channels' response. K-L transform, also known as 
PCA (Principle Component Analysis), can be used to characterize the independence of different 
channels since it de-correlated components of X{n) as completely as possible [5,6]. This is achieved 
by analyzing the eigenvalues of X(n)'s covariance matrix C. More evenly distributed eigenvalues 
would suggest more independence between the channels. The covariance matrix is defined as: 

N 



C = E 



(X-X)\X-X)'\x = ^Y^X{n) (4) 



where ^[.] denotes the expectation. Assume that matrix C's eigenvalue vector is K = W,A^,"'X^^^Y 

k In 

and \ >X2>"'> >^84 . PCA accumulative curve is defined as f{k)=^^X. ^ A. . The ideal / {k) is 

i=\ I i=\ 

a 45° straight line which goes through the origin point because all X^ are equal. In reality, however, 
/ {k) is a convex curve because of the dependency between the channels. 

Table 3 shows the five experiment cases that were carried out. Case 1 and case 2 would demonstrate 
independence of channels in float and mode 0. Case 3 to 5 are used to illustrate the utility of 
accumulative curve method since in cases 3 and 4, common X-ray flux noise and 50 Hz power noise 
are included respectively, which will affect all channels simultaneously and will impair independence 
quality considerably. These tests were also carried out using a single detector unit with the setup 
described in Section 2.2. 
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Table 3. Experiment cases of independence test. 



Case 


Test environment 


Frames N 


integration time 

T(us) 


Readout mode 


X-ray 


1 


Testl 


476 


500 


float 


Off 


2 


Testl 


456 


500 


modeO 


Off 


3 


Testl 


468 


500 


modeO 


90 kv,80 mA 


4 


Test2, 50 Hz power noise 


512 


500 


modeO 


Off 


5 


Test2, 50 Hz noise removed 


472 


500 


modeO 


Off 



2.3.3. Linearity and Sensitivity 

In the experiments we checked two situations. In situation 1, tube voltage which was denoted as 
U^^j^^ and integration time which was denoted as T were fixed, whereas tube current which was 
denoted as /^^^^ varied. In situation 2, t/^^^^ and /^^^^ were fixed but T varied. The objective is to 
observe the Unearity property with respect to /^^^^ and T, Tables 4 and 5 show the experimental 
conditions for each situation. 



Table 4. Situation 1, r= 500 us, Utube= 90 kV. 



Itube 


10 mA 


20 mA 


32 mA 


40 mA 


50 mA 


63 mA 


80 mA 


100 mA 


125 mA 


frames N 


72 


72 


72 


72 


72 


72 


72 


72 


72 



Table 5. Situation 2, hube = 80 mA, U,ube = 90 kV. 



T 


300 ms 


400 us 


500 MS 


600 us 


700 ms 


800 us 


frames N 


72 


72 


72 


72 


72 


72 



Collected data were denoted as xj{n), n = 0,l,L ,N , i is the index of channel i = l,2,L ,384 and 
j is the index of T or /^^^ . The linear model is showed as follows: 

y = /3.x + £^ / = 1,2,L ,384 

n=0 

where y is the response variable, x stands for T or /^^^^ . Parameters of J3. and £. is got through 

Y N-l 

linear regression by using least square method [7]. The mean value x.= — Vx/(/7) is used as 

regressor value. 

Sensitivity can be defined as: 

AR. 

^i^^tube)= .rj^'j (6) 

where Ai?. is channel's response increase, Ar is integration time increase and A/^^^^ is tube current 

increase. The physical meaning of sensitivity is the response increase of 1 mA increase of tube current 
and lus increase of integration time. Actually in experiment, sensitivity can be get through two 



Sensors 2013, 13 



5174 



methods: Method 1: s'j.=— and Method 2: s'j and their results are expected to be equal. 

^ ^tube 

Actually, Medhod 1 and Method 2 are two ways to get channel / 's sensitivity, however, their results 
may be slightly different due to inevitable experiment error. 

2.3.4. Uniformity and Uniformity Correction 

Due to the manufacturing process, the sensitivity of each channel in the detector unit may be 
different, that is s. ^ Sj y , which means for the same X-ray dose, response of different channel 
may be different. In order to evaluate the uniformity of sensitivity, we define the uniformity parameter 
Uni based on the responses. 

Stepl: acquire a fi*ame of data, denoted as x(i)J = 1,2,L ,384 ; 

Step2: x(i) is normalized to = x(/)/ max x(0 , where maxx(/)= max{x(0}; 

l</<384 

Step3: calculate Uni by the following procedure: 

383 ^ 383 

f«'W) = ^58^gWO-?f.J=5ijZi<0 (7) 

When the response is ideally uniform, s. are identical, so Uni = 0 . For a chessboard image, 

however, we should get Uni = 0.5 . Because of the non-uniform sensitivity, raw data collected fi*om 
the detector unit should be corrected before reconstruction. Poor uniformity will introduce ring 
artifacts into reconstructed image. One possible solution is to tune different sensitivity s to an identical 
value. The correction procedure can be written as: 

= ^/ag(^ , ^ , • • •, ^) • (7 - s) /g) 

^2 *^384 

where Y is raw response vector, e is the vector of £. in linear model (7), diag(L ) stands for a 
diagonal matrix and R^is the corrected data. Actually this procedure makes a simple linear projection, 
which tunes every channel's sensitivity to the maximum value among them, and therefore linear 
property of response is kept. The benefit of this transformation is that it makes full use of dynamic 
range. The drawback of the projection is that noise may be amplified as well. Considering the 
drawback of this procedure, a similar method may be used, which is: 

.S- S- s 



ct ft e ^ ^ 



^1 ^2 ^384 

where s^.^ is the minimum value of s. . The fact is that, however, the raw data is largely uniform and 
only a few channels are abnormal. Considering this fact, a more reasonable method is: 

R=diag{—,— 

(10) 

1 384 ^ ^ 
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This method reconciled drawbacks of Equations (11) and (12) and may yield a probabilistic 
good result. 

The former uniformity correction method is based on knowing each channel's sensitivity s. and s. 
varies with tube voltage t/^^^^ . Physically, t/^^^^ determines the X-ray's energy spectrum and thus, the 
relationship between s. and U^^j^^ is non-linear. In order to gain an adaptive correction method, we 
should test the s.-U^^^^ line through experiment and arbitrary s.(U^^^J can be obtained through 
interpolation. 

In the experiment, we used Method 7, that is s\- 15\ IT to test all channels' sensitivity s\ . Based 
on the results obtained, we corrected some randomly selected frames of data and a comparison of 
uniformity Uni between raw data and corrected data was shown. Finally, we showed a classical 
s. - U^^j^^ curve and the interpolated result. In order to show the s. - U^^^^ line, we repeated the 
procedure of testing sensitivity s'j. in Section 2.3.3. We tested the response of 10, 32, 50 and 80 mA at 
each tube voltage. The tested tube voltage is 60, 70, 80, 90, 100, 110 and 120 kV. The interpolated 
data is at 115 kV. 

3. Results and Discussion 

3.1. Noise Level 

Figure 6 shows the noise level of each pixel by calculating a and ju in Equation (2). During the 
calculation, was set to 60 for integration time 500 us, which was obtained from the manufacturer. 
Each point in the figure stands for a channel. Squares are channels tested in float mode, whereas the 
dots are channels tested in mode 0. These results show that the noise level in mode 0 and float mode 
meet the requirement since maximum standard variance maxcj is under 70 and most of dots and 
squares concentrate near (j = 55. In addition, mode O's dots are dispersed wider in axis (J compared 
with that of float mode because of the introduction of PDA noise . Furthermore, an obvious offset 
effect is introduced in axis jU . Obviously, the dominating noise source is . 

Figure 6. Noise level of float mode and modeO. 
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3.2. Channel Independence 

Figure 7 shows the PCA accumulative curves for single unit tests {test 1) and CT detector test 
(test 2) respectively. Generally, upper accumulative curve's independence is poorer compared with the 
lower one. For example, the line of testl, float mode is lower than that of testl, mode 0 since the 
introduction of n^ to the latter, impairing the independence quality. Another phenomenon illustrates 
this property clearly is that when we removed the 50 Hz power noise from the raw data in test2, its 
accumulative line dropped dramatically, which is showed in Figure 7. The rationale of this 
phenomenon is that power noise acts as a common noise source which affects each channel 
simultaneously, in another word, channels are highly dependent. In order to further explain the utility 
of this method, we introduced another common noise source — X-ray tube noise /2^^^^ . Usually n^^^^ is 
significantly more conspicuous than that of n^+n^^ n^ so that we can ignore the effects of the later 
when X-ray tube is turned on. In Figure 7, the n^^j^^ curve is also rather high, accord with our analysis 
and expectation. 



Figure 7. PCA accumulative curve. 




Table 6 shows parameters related to the tests. Samples is the number of frame . Noise mean is 

N 

mean value of all channels noise level ( O"). First component ratio \^ \l ^Y^^i and it is important since 

it reveals the independence quality directly, the smaller the better. Test environment is some important 
test information. 

Table 6. Independent test parameters. 



Samples Noise mean First component ratio Test environment 



Testl, float mode 


476 


57.0 


15.85% 


0 kV, 0 mA, 500 us, float 


Testl, modeO 


456 


57.8 


15.22% 


0 kV, 0 mA, 500 us, modeO 


Testl, quantum noise 


468 


7,831.3 


65.72% 


90 kV,80 mA, 500 us, modeO 


Test2, 50 Hz noise 


512 


80.8 


68.42% 


0 kV, 0 mA, 500 us, modeO 


Test2, 50 Hz removed 


472 


60.7 


1.37% 


0 kV, 0 mA, 500 us, modeO 
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3.3. Linearity 

In order to check the linearity of response of the detector unit, we randomly choose three channels 
and test its R-T and R-If^h^ line, T is integration time, I^^^^ is tube current and R is the response. 
Figure 8 (a,b) shows R - T and R - /^^^^ line respectively, which exhibit rather linear responses. 
Figure 8(c) shows channel 73 's sensitivity that was got using the 2 different methods described in 
Section 2.3.3. The blue solid line's slope stands for Sj and the red line stands for Sj. . Just as we have 
expected, they are almost the same. Figure 8(d) shows relative difference between Sj and Sj. of each 
channel, the maximum difference is 1.67%. 

Figure 8. (a) Linearity property of R - T , X-ray tube parameter is 90 kV, 80 mA; 
(b) Linearity property of R -I^^i,^ , integration time T= 500 us; (c) Sensitivity property, red 

line's slope is Sj. and blue line's slope is Sj ; (d) Difference between Sj and Sj. 




Tfbiue) or I K SOO/1 2S(r^d) No. channet 



3.4. Uniformity and Correction 

Firstly, we calculate the sensitivity s^ of each channel (since Sj is largely the same as s^ , we plot 
Sj here). Figure 9(a) shows the map of s^ . From this picture, we can find dramatic turbulence of 
sensitivity, which reveals poor uniformity character of PDA assembly. It is obliged to correct raw data 
before reconstruction. 
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Figure 9. (a) Sensitivity s^ . (b) Raw data and corrected results. The blue line is a frame 

data of 90 kV, 32 mA, green line is 90 kV, 80 mA, red lines are corrected results, 
(c) Corrected result of a randomly selected channel, obviously, this correction procedure 
keeps linearity, (d) Uniformity Uni, which is calculated by Equation (10). Red line is the 
corrected result, blue line is raw data's uniformity. 
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We try to correct some randomly selected frames of data by using the mean method of 
Equation (10), Figure 9(b) shows two results, which are two frames of data sampled at 90 kV, 32 mA 
and 90 kV, 80 mA respectively. Apparently, corrected data is far smoother than that of raw data. 

Actually, this uniformity correction procedure is a simple linear transformation so that response's 
linearity is kept. Figure 9(c) demonstrated this analysis. We linked the corrected result of a channel 
directly with line segment and the whole line is just like a complete line segment. Table 7 shows 
channel 371 's raw data and corrected data. 



Table 7. Channel 371 s raw data and corrected data, 90 kV, 500 us. 



I tube (mA) 


0 


10 


20 


32 


40 


50 


63 


80 


100 


125 


Raw data 10*' 
Corrected 10*' 


0.016 
-0.024 


0.414 
0.370 


0.808 
0.753 


1.302 
1.250 


1.626 
1.571 


2.034 
1.980 


2.577 
2.515 


3.262 
3.207 


4.066 
3.996 


5.053 
4.975 



Finally, we show the a contrast of uniformity Uni , which was defined in Equation (7). For almost 
all response, the corrected result has better uniformity. For the zero input, however, the uniformity is 
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undermined by the correction procedure, this is because raw data of zero input's uniformity is quite 
well, just as we have analyzed in noise level experiments. 

Based on this fact, one solution is that if raw response is too small, near the zero-input response, for 
example, correction should not be performed. This strategy can avoid the obvious error introduced by 
correction procedure. 

Figure 10 shows one classical s^-U^^j^^ line, which belongs to a randomly selected channel 313. 
Every channel has largely the same shape. We test sensitivity s^ at X-ray tube voltage of 60 kV to 
120 kV, covering clinical CT voltage range. Sensitivity at other voltage can be got through line or 
spline interpolation. The spline interpolated value at voltage 115 kV is 142.23 and tested value is 
141.92, almost identical. 

Figure 10. s^ -U^^^^ line of channel 313. 




60 70 80 90 100 110 120 

kv 



3.5. Reconstructed Images 

In whole CT detector test environment, we scanned some phantoms, such as an apple, a head 
phantom and a line pair phantom. The reconstruction method is fan beam FBP algorithm, without data 
correction. The result is showed in Figure 1 1 . The general inner parts of these phantoms are clearly 
showed in the image except for some detailed information. Since lack of uniformity correction, 
circle-like artifacts is rather obvious, especially in apple image and line-pair image. Stripe artifact is 
conspicuous in head phantom image. Further analysis of image quality is beyond the scope of this 
article, we will demonstrate more analysis and result in the future articles. 
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Figure 11. Reconstructed images of phantom, (left) slice of an apple; (middle) head 
phantom; (right) line pair phantom. 




4. Conclusions 

In this article, we described a modular detector unit which can be easily tiled together to form a CT 
detector. Experiments were conducted to characterize the detector performance both at the single unit 
level and system level. The noise level, linearity and uniformity of the proposed detector unit were 
reported and results suggest that it is satisfactory for CT imaging. Initial imaging studies were also 
presented which demonstrated the potential application of the detector unit in actual CT scanners. 
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